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1 Introduction 



The study of exactly solvable models has been an important activity in equilibrium statistical 
physics and it is equally desirable to find such models in the field of non-equilibrium phe- 
nomena. In this context it has proved very useful to exploit the formal analogy between the 
master equation for classical stochastic systems and the quantum-mechanical Schrodinger 
equation. For example, by relating the one-dimensional kinetic Ising model to a quantum 
spin chain |], |], §, |J it becomes obvious why the system is solvable with Glauber's || choice 
of the rates. Similarly, the hopping of classical particles with hard-core repulsion on a lattice 
can be formulated as a spin problem, namely the ferromagnetic Heisenberg model |], |7], and 
in one dimension this is again a solvable model. 

In recent work, this approach has been extended to more general situations. They in- 
clude diffusion with a preferred direction as well as reactions between particles or coupling 



to external reservoirs j8L 0, [10, 11, 12, 13, 14]. The spin of the resulting quantum models is 



s + 1 if s types of particles are involved. Fermionic or bosonic representations have also been 
used, for example for models of self-organized criticality [|T3 , |T6|j . In general, the inclusion 
of more processes makes the time evolution operator more complicated. On the other hand, 
by tuning the various rates one may also achieve simplifications which make the problem 
integrable in some cases. 

One of the systems where this happens is the coagulation model in which particles hop 
on a one-dimensional lattice and, in addition, can merge and separate again with certain 
probabilities. This model was introduced and treated in a continuum approximation by Do- 



ering et.al. |jTT|, [TBI, O, GDI]. It was found later that it is exactly solvable on a lattice if the 



rates for hopping and for coagulation are equal The time evolution operator can then 



be brought into the form of the Hamiltonian for a spin one-half XY chain with a Z field. It 
is the dual of the operator for the Glauber model mentioned above and can be diagonalized 
in terms of fermions. The spectrum has a gap determined by the decoagulation rate. If 
this rate vanishes, one finds an algebraic decay of the concentration with exponent minus 
one-half. The model can then be mapped onto a pair annihilation model [^T] and there 
is also an interesting physical realization in terms of excitons in the quasi one-dimensional 
substance TMMC J22|. 

In treating this model, it turned out that physical quantities like the density are most 
easily calculated from the probabilities to find empty intervals of arbitrary length ("holes"). 
Their equations of motion form closed sets and can be solved for a ring as well as for a chain. 



They have been used also in related models [23]. An intriguing aspect is that these closed 



sets even exist if the creation of particles is allowed ||24j| . The relaxational spectrum of the 
hole functions then has the form of the Wannier-Stark ladder found for lattice electrons in 



a homogeneous electric field fl2"5] , |2"E[ . However, no fermionic solution is possible in this case 
and the question of complete integrability has remained open so far. 

In the present article we study the properties of the coagulation model for the case of a 
preferred direction. The geometry will be the open chain where an interplay between the 
boundaries and the spatial asymmetry occurs. Our aim is to pursue the question of integra- 
bility further, and also to find the basic physical features of the model. We focus on the case 
where all processes have the same asymmetry and the rates for hopping and coagulation are 
chosen equal. Then the model contains two free parameters. We show that, as for symmetric 



1 



rates, the spin operator describing the time evolution can be brought into a form where it 
is quadratic in fermions. Thus one finds complete integrability also in this case. However, 
in contrast to simple hopping or symmetric coagulation the problem has no quantum group 
symmetry. This can be related to particular boundary terms which correspond to fields 
in the magnetic language and lead to special selection rules. There is, however, a discrete 
S^-symmetry which shows up in the multiplicity of the relaxational modes. 

We also consider the hole equations, both in the one-hole sector and in general. Following 
an approach by Bedeaux et.al. |27[] the different levels in the hierarchy can be decoupled and 
a complete solution achieved. By considering more general rates and processes, it is seen 
that this procedure only works in the free fermion case. This is also what the Reshetikhin 
criterion for integrability |28| gives. Thus the other cases where the one-hole equations can 
be solved, correspond to partially integrable systems. Such a situation occurs also in kinetic 
spin models P, 1211 and in reaction-diffusion models where for particular rates the equations 



for the density correlations form closed sets [30 



The physical properties of the model are also interesting. In the stationary state the 
density is inhomogeneous and determined by a competition between asymmetry and deco- 
agulation. If the latter is small, the density is finite only near one boundary. In the opposite 
case, it is constant in the bulk of the system with some additional boundary effects. In 
between, a transition takes place, where the profile becomes linear over the major part of 
the system. These different regimes are also reflected in the form of the relaxation spectrum. 
In general, one finds a gap but at the transition it vanishes and one has slow modes as in 
symmetric diffusion. The phenomenon may be viewed as another example of a boundary- 



induced transition JT3|, |3l], [32], j33f . 

We shall present the material as follows. In section 2 we introduce the model and the 
spin formalism used in its treatment. The time evolution operator is derived and brought 
into convenient forms in section 3. In section 4 its spectrum is found by diagonalization via 
fermions. Section 5 deals with the equations of motion in the one-hole sector and in section 6 
the n-hole problem is treated via cumulant-like functions. From these results, density profile 
and correlation functions in the steady state are obtained as special cases and discussed in 
section 7. Finally, section 8 contains a summary and concluding remarks. 



2 Model 

We consider classical particles on a chain of L sites, each of which can either be empty 
or occupied with at most one particle. The dynamics of the system consists of transitions 
which involve only the configurations of neighbouring sites. As discussed previously [p4 ], 



one can then distinguish altogether 12 different processes, namely hopping, coagulation and 
decoagulation as well as birth and death of single particles and of pairs. Each process is 
described by a rate constant and, in the absence of detailed balance, all these constants are 
independent. 

In the bulk of the paper we will only consider the first three processes which are shown 
below together with their rates: 
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The left-right asymmetry which we assume could result from external fields. We will use 
a spin language to describe the system. The state of site n is specified by a variable a n = ±1 
such that <j n = +1(— 1) if the site is occupied (empty). A configuration of the total system 
is denoted by a = {01, cr 2 , . . . cr L }. The probability P(a, t) to find configuration a at time t 
then obeys the master equation 

^P(^t) = "E P(^t)w(a^ a!) + £ P(a!,t)w(a! - a) (2.1) 

ct' a' 

where the two terms on the right represent the loss and the gain processes, respectively. 
The transition probabilities w(g_ — > a!) are sums of contributions from the various nearest 
neighbour processes. 

To formulate the problem in quantum-mechanical terms ffl, Q one uses a ket vector 
notation in the 2 L -dimensional state space, writing 

\P(t)} = J2P(a,t)\a). (2.2) 

The master equation then takes the form 





Of 



P(t)} = -H\P(t)) (2.3) 



from which the analogy with the Schrodinger equation is obvious. The time evolution opera- 
tor H will therefore be called Hamiltonian in the following. However, since \P{t)) is already 
the probability, expectation values differ from the quantum-mechanical ones. For a quantity 
A(a_) one has 

(A(t)) = J2A(a)P(a,t) = (0\A\P(t)) (2.4) 

with the bra vector 

(0| = £ (a\ . (2.5) 
The operator H is the sum of nearest-neighbour terms 

L-l 

H = } j H n ^ n+ i (2.6) 

n=l 
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and the (4 x 4) matrix H njn+ i can be written down easily |[24|| . In a basis of states ( | ), | — 

+), | H — ), | + +)) the result is, including birth processes (rates 6r,&l) and pair creation 
(rate c) for further reference 



H 



n,n+l 



( c + b L + b R 

-b R 

-b L 
V -c 





— 0£ 

-d L 



^ 

an + d R -c L 
—djR c L + c R ) 



(2.7) 



Due to probability conservation, the elements in the columns sum to zero. H n ^ n+ i and 
therefore the total operator H is non-hermitean. 

In the next sections we will consider the case of no birth and pair creation (b R = b^ = 
c = 0), equal rates for hopping and coagulation (c R = a R , Cl = cll) and equal asymmetry 
for all processes. Up to a common constant which we set equal one, the rates can then be 
written ^ 

a R = -, a L = q, d R = — , d L = Aq . (2.8) 

q q 

Thus the parameter q as usual describes the asymmetry, while A gives the ratio of decoag- 
ulation to coagulation rates and therefore governs the (average) density. 



3 Time evolution operator 

Starting e.g. from Eq. ( |2.7| ), the Hamiltonian H can be expressed in terms of Pauli-matrices 
(a=x, y,z). It is the sum of diffusion and coagulation contributions 

H = Hd + He , (3.1) 

where Hp is given by 

H d = -EW^i + ^^i + fe + f 1 )^! (3-2) 
4 1=1 L 

+(9-r 1 )«-< + i)-(g + <T 1 )}- 

The standard ?7 g [S'f/(2)]-symmetric form of Hd is then obtained by a position-dependent 



rescaling of a + and o R O, |TJ]. Physically, this corresponds to factorizing out of P(a,t) 

1 /2 

the quantity P (a), where Po(s) is the stationary distribution in the case of simple hopping. 
The coagulation part reads 



H c = ~\ £{(<TV + Ago-+)(l + a* +1 ) + (1 + ofi(qa^ + Ag'VtO 
4 i=i L 

- (g + g- 1 )^ - A) ^ - (g + g" 1 )^ + ^ +1 ) (3.3) 
+ A(g-g- 1 )(a l 2 -af +1 ) - (g + g- x )(l + A)} . 

It is more complicated than Hd because of the single spin terms af, a^, a\ which appear in 
it. One should note, however, that the contribution (a* — of +1 ) is effectively only a boundary 
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term. For q — 1 it vanishes altogether and H reduces to the expression found previously 
[H]] . In that case, H was brought into a simpler form by two local transformations in spin 
space: A hyperbolic rotation around the z-axis with generator 

V x = I] e A< ; e- 2X = VA (3.4) 

n 

which rescales the cr ± operators, and a real rotation around the y-axis with generator 

V 2 = l[e-^; co S (6) = -=L= (3.5) 

n V 1 + A 

which forms linear contributions of o x and cr z . 

It turns out that the same transformation can be used also for q ^ 1. The result for the 
new operator H' = VyVxHY^V^ then is 



D L ~ l ( I - a 2 



= -7 E + l~ lcr i a i+i + i 

Z 8=1 L 



l + q 

2 



r(f^ + q*Ui) ~ (7 + 7- 1 ) (3.6) 



1-g 2 

1 +q 2 



1~1 1 



1 



where 7 = a/1 + A. and Z) = |(g + q , ~ 1 )7- 

This operator has two non-hermitean parts: the third term, which is a kind of chiral 
contribution from the bulk, and the last term which again acts only at the two boundaries. 
Both vanish in the symmetric case (q = 1) and H' then reduces to the dual of the Glauber 
model operator. 

The bulk terms in H' are all quadratic in fermion operators after a Jordan- Wigner trans- 
formation. There are, however, also boundary terms in a x , a y and a z . The last one is 
unproblematic, but a x and a y are linear in the fermions. There is, however, a well-known 
way around this difficulty: One extends the chain by one site at each end and uses the addi- 
tional spin operators, which are constants of the motion, to make boundary terms quadratic. 
Thus the operator H' is solvable via fermions. 

Actually, one can bring it into a slightly simpler form by performing a rotation by 180° 



V 3 = J] (3.7) 

n 

and a further transformation with the generator 



to obtain 

L-l 



H 



\l E((^ + TVf^ + i + i + (3-9) 

z i=i 1 

+ qat +1 ) ~ \{q - - 7 _1 ) W " <&i) ~ (V + V' 1 )} 
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where (77 + r]^ 1 ) = \(q + q~ l )(^ + 7 _1 ). 

There is no o v boundary term in this expression so that the only non-hermitean part is 
the chiral bulk term. This operator will be diagonalized in the next section. 

We mention that H' and H" have a simple symmetry, namely they are invariant under 
the tranformation q 1/q combined with spatial reflection i — > L — i + 1. 



4 Analysis of the spectrum 



In order to diagonalize the Hamiltonian ( |3.9|) in terms of fermions, we write H" = H{f + H" 
where 

tff = ?K-^) (4-1) 



is the linear boundary term and 



^(?-ff- 1 )(7-7" 1 )- (4-2) 



As outlined in the last section, one appends one site at each end of the chain and extends 
H" by two Pauli matrices ctq and Then 

H'l = %lat-alal^) (4.3) 



is bilinear in fermions so that standard diagonalization methods |35) can be applied^. The 
eigenvectors of the extended Hamiltonian H" = Hq + H" include those of the original chain 
H" as follows: Since the appended matrices are constants of the motion, the spectrum of 

H" decomposes into four sectors (++, H — , — h, ) corresponding to the eigenvalues of <Tq 

and This implies that all eigenvectors of H" have the form 



® m ® , (4-4) 





where is a vector with 2 L components. Since in each sector the appended matrices can 
be replaced by their eigenvalues, the spectrum of the original chain appears in the (++) 
sector and the vectors are just the eigenvectors of H" . Therefore we can determine the 
eigenvectors of the original problem by projection onto the (++) sector. One should mention 
that a very similar situation occurs in calculations using the corner transfer matrix of the 
two-dimensional Ising model |36|, [37|]. In that case the sectors arise from the fixed inner and 
outer spins in the corner geometry. 

Defining anticommuting Clifford operators Tj' v = (Yl{Zo a i) a j' V one can rewrite H" as a 
bilinear expression in rf and rj. Then a second linear transformation = | Y^=\{.4>%^ ± 
io'Ly'') yields 

H" = (E A * fl K) - K (4-5) 



3 In the same way one could also add c^-matriccs or any linear combination of a x and a y . 
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which is a diagonal expression in fermionic creation and annihilation operators a k and a k . 
The one-particle energies A& and the eigenvectors 4> k = ((f>%j, ±z^ jl; . . . , <f>%,L, can 
be derived from the eigenvalue problem M^> 



± 



M 
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-D T 



V 



D 
B 

-A T 



A 
B + C 



A 

-A T 



B + C 
—A T 



A 
C 
D T 



-D 



(4.6) 



where A, B, C and D denote the 2x2 block matrices 



A= ?), B 

\ 7] + 77 iq 



Y-l 



,-1 







c 




D 




The eigenvalues of M are real and occur in pairs with different signs. It should be emphasized 
at this point that although a k 

{a k , aj } = 5ku we have 7^ {%) because of the nonhermiticity of H" 

Assuming that the wave function 0^ consists of exponential contributions e~^, the bulk 



and a k obey the usual fermionic commutation relations 



equations (|4.6|) yield the dispersion relation 



A 



7 (rj + 77 1 — qe 



(4.7) 



Including the boundary conditions which determine (, we obtain three types of fermionic 
one-particle states: 



• One trivial zero mode (Ao = 0) which is localized at the ends of the extended chain. 
This zero mode represents conserved quantities <Tq and o"£ +1 and causes two-fold de- 
generacies of all levels. 

• Two exponential modes with excitation energies A/ = Xjj = k. Their wave functions 
(pf and (f)jj decay exponentially with three different length scales corresponding to the 
values 111(7), m (? 2 7) an d ln(g 2 7 _1 ) for (. These modes reflect the influence of the 
boundary term H'{. 

• L — 1 oscillatory modes with A^ = 7(77 + T]~ l ) — 27 cos ^ where k = 1, . . . , L — 1. The 
corresponding wave functions show damped oscillations of the form <pk,j ~ g±i e ±l7rfc / L . 
This implies the presence of a fourth length scale with ( = Inq. The damping is a 
consequence of the nonhermiticity of the Hamiltonian for q ^ 1. 



One can show that the lowest level of H" is a two-fold degenerate state in the sectors (H — ) 
and ( — h) with a negative eigenvalue A = —k. The lowest level in the (++) sector with 
eigenvalue zero is already an excited state. The total number of excited modes in this sector 
has to be odd and one can distinguish: 



7 



• States with an even number of excited oscillatory modes. The levels are two-fold 
degenerate because either mode / or mode // can be occupied. This class contains the 
ground state with A = and the lowest gap is 

AA = 7(77 + --2). (4.8) 
V 

• States with an odd number of excited oscillatory modes. Here the modes / and II 
are either both occupied, which gives a contribution to A, or both empty, which 
gives (— k), i.e. we have two bands of excitations which are shifted symmetrically by 
±/t. Their lowest gaps are 

AA(k) = 7(7? + — - 2) , AA(-k) = 7 (- + - - 2) . (4.9) 

7g q 7 

The nature of these three branches (singlets and doublets) suggests a symmetry in the 
problem. This is indeed the case: From the operators for the two exponential modes one 
can construct the quantities 

a = 1 — (afaj + afjajj) + (afajj + af T aj) (4-10) 

r = 1 - ^(afaj - afja n ) 2 - ^-{afajj - a+aj). 

They satisfy the relations 

a 2 = 1, r 3 = 1, ara = r 2 (4.11) 

which characterize the generators of the permutation group S3 and commute with both H" 
and H" . Thus the levels can be chosen as irreducible representations of 5*3 which, according 
to group theory p8fl , consist of two types of singlets and one doublet. The operator a has a 



simple meaning, it interchanges the two degenerate fermion modes. 

The smallest gap is obtained from one of the expressions ( f4.9| ) which differ by a reflection 
q — > q~ x . It vanishes if 7 = q (q > 1) or 7 = q^ 1 (q < 1). One band of excitations is then 
massless and algebraic long-time behaviour is expected. At this point a transition takes 
place between a low-density phase (for 7 < q) and a high-density phase (for 7 > q). This 
will be seen in section 7 from the steady-state properties of the system. One should note, 
however, that the situation is different from a Pokrovskii-Talapov type of transition |39| 
since here the spectrum has a gap on both sides (7 > q) and always varies quadratically with 
the momentum k. 



5 Hole formalism 

Although the diagonalization in terms of fermions solves the model in principle, another 
complementary approach is simpler for practical purposes and gives additional insight. In 
this approach the states of the system are described in terms of strings of empty sites 
(holes) instead of spin configurations [IT]. In this way one can distinguish different sectors 



characterized by the number of holes. The equations for the probabilities couple a sector 
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only to the one with lower order. For the one-hole sector of the coagulation model (|2.7|) with 
free ends they were already derived in Ref. [2H]. Denoting by Q(x,y,t) the probability to 
find the sites x + 1, . . . , y empty and using the abbreviations 

Q-L,R = cil,r + d L)R - b LtR - c (5.1) 
1l,r = a L + a R + d LtR 
8 = b L + b R + c 

they read 

• for holes which do not touch the boundaries (0 < x < y < L): 

^-tt(x,y,t) = a R tt(x - l,y,t) + a L Vt(x + l,y,t) (5.2) 
at 

+ a R Q(x, y - 1, t) + a L Q(x, y + l,t) 
- (il + 1r + (y - x - 1)5) Sl(x, y, t) 

• for holes touching the left boundary (0 = x < y < L): 

^-M0,y,t) = a R n(Q,y-l,t) + a L n(Q,y + l,t) (5.3) 
at 

-( lL + (y-l)S)n(0,y,t) 

• for holes touching the right boundary (0 < x < y = L): 

^-Vt(x,L } t) = a R n(x-l,L,t)+a L n(x + l,L,t) (5.4) 
at 

-( lR + (L-x-l)S)Q(x,L,t) 

• for the hole extending over the whole chain (x = 0,y = L): 

—n(0,L,t) = -(L-l)5n(Q,L,t), (5.5) 
at 

where we formally define Q(x,x,t) = 1 according to [p4|j . This is the coupling to the 
lower sector (no holes), which here leads to an inhomogeneous system. The integrable case 
discussed in the preceeding sections corresponds to 

«l,h = 7V 1 , a L:R = q ± \ Ti.^^+TV 1 , <5 = . (5.6) 
Introducing rescaled probabilities Q(x, y, t) via 

Q(x,y,t) = -yy-*if**n(x,y,t) (5.7) 
and separating the time dependence exp(-At), Eqs. ( |5.2|) - ([5.5|) reduce to: 

('Y(q + q- 1 )('y + 'y- 1 )-A)to(x,y) = 7 (fi(ar - 1, y) + Q(x + 1, y) (5.8) 



+tt(x, y - 1) + Q(x, y + 1) 

(^ 1 +7 2 g-A)^(0,y) = 7(fi(0,j/-l) + f2(0,j/+l)) (5.9) 

(9 + 7 2 f 1 -A)n(x,L) = -y(p{x + l,L) + tt(x-l,L)) (5.10) 

AQ(0,L) = (5.11) 
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with the inhomogeneous boundary condition Cl(x, x) = q 2x . These equations can be solved 
using similar techniques as in Ref. pT| which rely mainly on the invariance of the bulk 



equation ( |5.8| ) under reflections x y and x <->• L — y. A trivial inhomogeneous solution is 
the empty lattice, for which fl(x,y) = ^y-~ x qy +x . The homogeneous set of solutions satisfies 
the condition Cl(x,x) = and reads 

Qofay) = h(L-x)g(y)-h(L-y)g(x) (5.12) 

Ao = 

(n , . . nkx , , . . nky . . 

<Pfc {x,yj = sm — — /i(L — y) — sm — — h{L — x) (5.13) 

A[ L) = 7 ( 7 g + 7 - 1 g- 1 -2cOS^) 

*fe = sm—j—g(y)-sm——g(x) (5.14) 

Af } = 7 (7^ 1 + 7^-2cos^) 

< =Pfc,z( a; ? Z/) = sm — j— sm — sm — sm — j^- (5.15) 

A k>l = 7 ((g + g- 1 )( 7 + 7 - 1 ) -2(C0S^ + C0S^ 

where l<k<l<L — 1 and the functions g and /i are 

g(x) = M! - , hM = - w«>- . (5.i6) 

L Ft 

We can now compare the eigenvalues Ao,A fc ' and Aki with the spectrum computed in the 



last section. As in the symmetric case pl| , A gives a two-fold degenerate ground state 
while Afr' R and A^ are seen to be one- and two-particle excitations. However, in contrast 
to the symmetric case we have A£ ^ Af since these levels are shifted by ±«. Therefore the 
boundary effects of H'l do not change the selection rules, they only modify the mass gaps of 
the one-particle modes. 



6 Probabilities for n holes 



For a complete treatment of the system one has to study configurations with an arbitrary 
number of holes. This was already discussed in Ref. |24| and we will now carry out such a 



study. The probability to find the strings of sites [x\ + 1, • • • , yi), (x 2 + 1, • • • , y^)-, • • • , (x n + 
1, • • • , y n ) with < X\ < ■ ■ ■ < y n < L empty at time t is denoted by Q n (xi,yi, • • • , x n , y n , t) 
and will also be called n-hole function. If no hole extends up to the boundaries, the time 
evolution of Q n takes the form : 



d 
dt 



^ n (x 1 ,y 1 , 



*£ni Vr, 



r 

t) = -(jL+lR+iyi-Xi-tfojttnix^yi, 



y n ,t) 
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+a R n n (- ■ ■ , Xi - 1, yi, ■ ■ ■ , t) + a L n n (- ■ ■ , Xi + 1, y h ■ ■ ■ , t) 
+a R Vt n (- ■ • , yi - 1, x i+1 , •••,*) + a L n„(- • • , y { + 1, ••-,£) 



(6.1) 



with the parameters defined in ( |5 . 1| ) . The i th term in the sum on the right hand side can 
be identified as time evolution of a one-hole function ( |5.2|) . In the following we choose the 
reaction rates according to (2.8) and rescale the fl n as in the previous section : 

Cln(xi,yi, • • • , x n , y n , t) = 7 Eiid'i-*0 gEiifo+w) ^(a^, yi , • ■ • , x n , y n , t) . 



(6.2) 



Then the rescaled time evolution equation turns out to be invariant under an arbitrary 
permutation of the variables (x\, ■ ■ ■ , y n ). This symmetry can be taken into account by 
introducing new coordinates (zi, z 2 , • ■ ■ , z 2n ~i, z 2n ) '■= {x\, yi, • • • , x n , y n ). The complete sys- 
tem of equations then reads: 

• if no holes touch the boundaries (0 < z±, z 2n < L): 



d_^ 
dt 



tt n (z u - ■ ■ ,z 2n ,t) = -717(7 + 7 l )(q + q l )Q, n {z x , ■ ■ • , z 2n , t) 



(6.3) 



2 ii 



+ 7E 



a 



■ > z i ~ 1) z i+li ' ' " > t) + ^n(- • ■ , Zi + 1, 2^+1, • • ■ , t) 



if the leftmost hole touches the boundary (0 = z%, z 2n < L): 

j t & n (o,z 2 ,---,z 2n ,t) = -((n- 1)7(7 + i~ 1 )(q + q~ 1 ) + i 2 q + q~ 1 )& n (zi, 



2n 



+ 7E 



i=2 



^n(' ' ' > z i ~ 1) ' ' ' ) t) + ^n(' ' ' j z i + lj z i+\i ' ' ' j t) 



if the rightmost hole touches the boundary (0 < zi, z 2n = L): 



d^ 



dt 



^n( z ir- -,Z2n-i,L, t) = -((n- 1)7(7 + 7 + ? 1 )+7 2 <? + <? ^n^l, 



2n-l r 



+ 7 E 



i=l 



( ' ' ' , Zi 1 , Z> 



j Zi ±, , 



£) + f2 n (- • • , 2j + 1, Zi+i, ■ ■ ■ , t) 



if the leftmost and the rightmost holes touch the boundaries (0 = Zi, z 2n = L): 

■^-Q n (0,Z2, ■■■,Z 2n -l,L,t) = -(n- 1)7(7 + 'J^iq + q^&n^ir ■ ■ , z 2n,t) 

dt 



2n-l 



+7E 



i=2 



\ ' ' 1 Zi 1 , Z> 



5 z i+ll ' ' ' i 



t) + a 



nV ' ' ) z i + 1) z i+li 



z 2n, t) 

(6.4) 

(6.5) 
(6.6) 



However, this set of linear equations for the n-hole functions is not closed since on the right 
hand side terms with Zi = Zi+i appear. These terms have to be identified with (n — l)-hole 
functions according to 



&n{ z li ' ' ' j z 2n) — q 2Zl &n-l{ z l> • • • , z i-l, z i+2, " " ' , z 2n) if z i — Z 



i+1 



(6.7) 



and therefore cause a coupling to the lower sectors as mentioned before. This coupling 
constitutes the main difficulty in solving the equations. A similar situation occurs for the 
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spin correlation functions of the Glauber model, and in that case a method to solve the 
problem was developed by Bedeaux, Schuler and Oppenheim |27| . It can be applied in the 
present case, too. For that purpose we introduce new cumulant-like functions C n with the 
property that 



(1) C n obeys the differential equations (|6.3|) - (|6.61 ) 

(2) C n [Zi, • ■ ■ , Z2n, t) — if Zi — Zi+i 



Following Ref. [27], these functions are defined by 

c n (z ir --,z 2n ,t) = £(-i) |p| (-i) (fc - 1} (k-i)\n ni ( z (, 



z" 

5 z 2n 



X.*) 



(6.8) 



Xfl n2 (^2ni+l> ' ' ' ' Z \n x +2n 2 i t) ' ' ' ^n k { z 2n-2n k +li * ' ' > Z \ni ■ 



The sum on the right hand side of (|6.8| ) runs over partitions £ of the numbers Z\, • • • , z 2n into 



with sizes rij where 



k Subsets (zf, • ■ ■ , Z2 Ul ), (^2rii+l) ' ' ' 5 z 2n 1 +2n 2 )^ ' ' ' i [ Z 2n-2n k +ly 

j = 1, • • • , k. Within these subsets the z\ are ordered according to their magnitude. P is the 
permutation which transforms (z%, ■ ■ ■ , z 2n ) into (zf, ■ ■ ■ , ^n) an d \P\ denotes its (uniquely 
defined) sign. The first two functions C\ and Ci read explicitly: 



Cl{Zi, Z2, t) 
Ci{Z\i %2, Z3, Z4, t) 



Q l (z 1 ,z 2 ,t) (6.9) 

^2(^1, z 2i z Zi z i, t) — Ql{Zi, Z2, t)Qi(z 3 , Z4, t) 

-Qt(zi, z^t)^}^, z 3 ,t) + f2i(zi, z 3 ,t)&i(z 2 ,Z4,t) (6.10) 



The proof that the C n have the property (2) by which the sectors are decoupled, is completely 



analogous to the one given in |27j. Property (1) holds since each term on the right hand side 
of equation ( |6.8D solves the differential equation (|6.3|)-(|^o]). 

However, this is only true if the conditions 5 = and olr(i l = a^a R are satisfied so that 
the rescaled equations d6.3| )-( pT6]) are symmetric in the starting points x and end points y of 
the holes. In the general case, i.e. for arbitrary values of the parameters 5, aL,aR,a,L and 
an, this symmetry is lost and arbitrary products of ^-functions are no longer solutions of 
the generalisation of ( |6.3| )-( |6T6| ). Then the given construction of the C-functions fails. It 
may not be so surprising that the method of Bedeaux et.al., which was developed for the 
Glauber model, works also here only in the free fermion case. To get an idea whether the 
general case is integrable or not, one can apply Reshetikhin's criterion for integrability |p8fl: 



H n + H, 



n+l, 



n+l 



X n — x. 



n+l 



(6.11) 



to the Hamiltonian ( [2.7|) , where X n is an arbitrary 4x4 matrix. We found that it leads to 
the same restriction on the reaction rates, namely 5 = and a^aL = oilcir- It is known, 
however, that this criterion is not a necessary condition. 

Turning to the solution of the equations, one first observes that they are inhomogeneous 
only for n = 1. A particular solution C p is again the empty lattice: 



C[(z u z 2 ) 



^(22-21)^(22+21) 



Ci 



for n > 2 



(6.12) 



12 



For the solutions of the homogeneous system we write 

C* ( Zl , • • • , z 2n , t) = e~ At & n \z u • • • , z 2n ) (6.13) 

and are left with an eigenvalue problem as in the previous section. As there, one finds four 
types of solutions (c.f. Eqs. (fngfl -(glg)): 

< k2 ,-M n = E(-l) M n^^) (6.14) 

reS2n »=1 

2n irk- 

A = 727(7 + 7 _1 ) (<? + Q 1 ) - 2 7 X! cos 

i=l L 



= E (-i) M (^-^(i))II^ -rMo) ( 6 - 15 ) 

reS 2n V i=2 7 -' 



2n 7rA;- 

A = ((n - 1)7(7 + T l ){q + q' 1 ) + 1\ + q' 1 ) - 27 E cos V 



,=2 i 



CS,,^ = E (-1)' T| II Biny, r( ^(, T(2n) ) (6.16) 

2n-l ^ 

A = {{n - 1)7(7 + i' l ){q + q' 1 ) + i 2 q + q~ l ) - 27 E cos V 

i=i ^ 



/ 2n-l u \ 

<J£U^ = E ("1) |T| [KL ~ MD) II sin V W(*(*o)) ( 6 - 17 ) 

reS 2 n V i=l ^ 



2n-l 

A = (n - 1)7(7 + 7~ 1 )(<7 + 'T 1 ) - 2 7 E cos ~r 

i=2 ^ 

In these equations the ki are integers with 1 < k\ < ■ ■ ■ < k 2n < L — 1. The sum runs over 
all elements of the permutation group S^n- To understand these solutions we remark that 
each term in the sum of ( |6.14| - |6TTT|) solves equations ( |6.3[ - |6.6D up to boundary conditions. 
These are satified by taking the totaly antisymmetric combination. In the case n — 1 one 
recovers the results of section 5. 

Comparing the eigenvalues A with the spectrum given in section 4 we find that they 
correspond to the excitation of 2n, 2n — 1 and 2(n — 1) fermions for the cases fl6.14|) , (|6.15| )- 
( |6.16j ) and ( |6.17| ), respectively. Furthermore, one has as many solutions as n-hole functions. 
Therefore the system of eigenfunctions found above is complete and the problem is thereby 
fully solved. 



7 Steady state properties 

The simplest application of the previous results is the investigation of stationary properties. 
The system has two steady states, namely a trivial one (the empty lattice) and a nontrivial 
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one where particles are present. Considering the situation where at least one particle is 
present (i.e. fl stat (0,L) = 0) we have: 



n stat (x,y) = l-^ L ~y^q^y-^%(x,y) (7.1) 

Inserting $ from Eq. (|5.12| ) then gives the following exact expression for the concentration 
c(j) = 1 - ttstatU ~ 1, j) at site j: 

c(j) = ^{l 2L (V - 1) + {q 2 - 1)7 2 (?7)- 2J ) - q 2L ((7 2 - 1)^ + {q 2 ~ l)(qh)^) } (7.2) 
where 

k = 7 2 (7 2i + 7" 2L - q 2L - q~ 2L ) (7-3) 

The density profile c(j) is shown in Fig. 1 for 100 sites, a fixed asymmetry q = 1.2 towards 
the left and various decoagulation rates A. As can be seen, there are three types of curves: 
If A is small, particles are found only near the preferred (left) boundary, while for large A a 
plateau exists together with some additional boundary effects. At the border between these 
two cases, the concentration decays linearly in the bulk. Computing the leading order of 
c(j) for L — > oo gives the following explicit forms for the profile 

[ £sf* + ferC*)"* if A < q 2 -1 

c(j) = ^(qi-V + l-l) if A = g 2 -1 (7.4) 

[j^ + iq 2 - l)((g 7 )- 2j - ^) 2L - 2f ) if A > g 2 - 1 

From these expressions one sees that three different inverse length scales given by 41n(g), 
2 ln(g/7) and 21n(g7) appear which are related to those found in the fermion eigenf unctions 
of section 4. The mean concentration per site c in leading order for L — > oo is 



(7.5) 



A 
L 


if 


A < 


q 2 


- 1 


1 A 

2 1+A 


if 


A = 


q 2 


- 1 


A 

1+A 


if 


A > 


q 2 


- 1 



where 

A - (l + A-g 4 )(l + A-g- 2 ) 

(1 + A-g 2 )(l + A)(g 2 -g- 2 ) " ' 1 ' ' 

Thus we have to distinguish two different phases. In the low- density phase A < q 2 — 1 the 
asymmetric diffusion is strong enough to move the bulk particles to the boundaries where 
they can coagulate. Here one has a stationary state with a finite number of particles A > 1 
and therefore the mean concentration c is of order 1/L. In the high- density phase A > q 2 — 1 , 
the bulk concentration yA^- is finite and independent of the asymmetry q so that c is of order 
one. At the transition A = q 2 — 1 (i.e. q = 7 or q = 7" 1 ), one of the lengths diverges and 
one has a linear decay of the concentration in the bulk. This is also the point where the gap 
in the spectrum vanishes (c.f. section 4). 

This phase structure also shows up in the stationary two-point correlation function. 
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Denoting by rii and rij the occupation numbers at the positions i < j, its connected part is 
given by 



g(i,j) 



(niUj) - (rii) (nj) 

((l - m)(i - nj )) - {{l - m)) ((l - nj )) 

n° tat (i - i,i, 3 - 1,3) - n? at (i - M)nf*(j - U) 

7 2 ( i-i-i) ^ n; « ( . _ i, j - 1) - nf*(i - i, j - m? a \i,jf) 



(7.7) 



To derive the last formula one uses Eqs. ( |6.9| ), (|6.10 ) and the f act that for stationary states 
all functions C n for n > 2 vanish because their relaxational spectrum has no zero eigenvalues. 
Inserting the expression for the funtion Ql tat (i,j) one has: 



g(i,j) 



L+j-i-3 qL-j-i+1 



( h(L - i + 1) - (7^) fc(L - i)) ((t?) ^(j - 1) - $0')) 

- (( 7 g) /i(L - j + 1) - - j)) (^(i " 1) " (T^ 1 ) 5 (i)) 
-c(z)c(j) (7.8) 

In the symmetric case (q = 1) the system has a ground state of product form. Therefore the 
particle concentration A/(l + A) equals the mean field value and the connected correlation 
function is zero. If there is an asymmetry to the left (q > 1), one finds the following properties 



• In the low-density phase A < q 2 — 1 the only region where correlations are present 
is the left boundary. Here the correlation function is negative since the dominating 
coagulation process reduces the probability that particles meet each other. 

• In the high-density phase A > q 2 — 1 there are correlations at both boundaries. As in 
the previous case they are negative at the left boundary and vanish in the bulk. At 
the right boundary where the particle concentration is reduced, decoagulation is the 
dominating process which increases the probability of neighbouring particles. 

• At the phase transition A = q 2 — 1 there is a linear decay of the correlation function 
for large distances. For large L and i/L, j/L fixed one has 

g(i/L,j/L) = ( T ^ z ) 2 (i/L)(l -j/L) (7.9) 

For short distances the correlations are nontrivial only near the left boundary where 
they have exponential form. For large L and i, j fixed 

9(i,j) = "(t^a) 2 (! + ^ _44+2 ) ^" 4J+2 ( 7 - 10 ) 

In this case the remaining length scale 1/4 111(7) enters and describes the boundary 
effect. 
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8 Concluding remarks 



We have studied the combined effect of a preferred direction and of open boundaries in a 
one-dimensional coagulation model. By formulating it as a quantum spin chain, the complete 
integrability was shown. The relaxation spectrum, although obtained from a non-hermitean 
operator, turned out to be real, with a gap determined by the asymmetry as well as by the 
decoagulation. An underlying S^-symmetry of the Hamiltonian was found and the relax- 
ational modes seen in the hole functions could be understood from the general result. 

The hierarchy of equations for these functions was considered and decoupled by a method 
introduced previously for the kinetic Ising model. In this way a complete solution was 
achieved. The decoupling approach lead to the same restrictions on the rates as Reshetikhin's 
integrability criterion. There are other solvable cases where a relation to the XXZ Hamilto- 
nian is used [|30|] . Whether further integrable situations exist, remains open. 

As to the physical properties, it was found that the model has two phases with low and 
high average particle density, respectively. The corresponding density profiles and correla- 
tions were calculated. They contain the various lengths which arise from the interplay of 
decoagulation and asymmetry. The approach to the steady state also turns out to be inter- 
esting. Calculations show that at the transition point an initially full system first develops 
a plateau at the mean field value of the density before a slow relaxation to the linear density 
profile sets in. Only at the boundary one finds a simple algebraic decay with a t" 1 / 2 law. 

A system with simple open ends, as considered here, is not the only possible case. More 
generally, one can supply and withdraw particles at the ends. Recent work has shown that, 
for pure hopping on the chain, this leads to interesting physical phenomena (boundary in- 
duced transitions) and to new mathematical features (matrix-product states) |]40| , |4l| , [14|, 
Since these processes are described by single spin operators at the ends of the chain, they 
can also be included in the present model. This is currently under investigation. Another 
possible direction would be to describe the coagulation process in more detail. This, how- 
ever, would lead to models with more than two states per site which are more difficult to 
solve. 
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Figure Captions: 

Fig. 1: Particle concentration c(i) for q = 1.2 and various values of A on a chain of 100 
sites. 
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Figure 1 



